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Abstract More than 80 giant planets arc known by 
mass and radius. Their interior structure in terms of 
core mass, number of layers, and composition how- 
ever is still poorly known. An overview is presented 
about the core mass M core and envelope mass of met- 
als Mz in Jupiter as predicted by various equations 
of state. It is argued that the uncertainty about the 
true H/He EOS in a pressure regime where the gravi- 
tational moments Ji and J4 are most sensitive, i.e. be- 
tween 0.5 and 4 Mbar, is in part responsible for the 
broad range M core = - 18 M®, Mz = - 38 M m , and 
M core + Mz = 14 — 38 currently offered for Jupiter. 
We then compare the Jupiter models obtained when we 
only match J2 with the range of solutions for the exo- 
planet GJ 436b, when we match an assumed tidal Love 
number fc 2 value. 
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1 Introduction 

Jupiter is the best studied giant planet. Its oblate 
shape and the corresponding elliptic gravity field de- 
formation have been measured within an accuracy of 
10~ 4 . The uncertainty of its mass reflects the un- 
certainty of the gravitational constant G, and even 
its atmospheric He abundanc e is known within 3% 
(see iGuillot and Gautierl 120071 for an overview) . How- 
ever, despite several spacecraft missions to Jupiter and 
observational data sampling over several centuries, we 
even do not know whether Jupiter possesses a core or 



Nadine Nettelmann 

University of California, Santa Cruz, Dept. of Astronomy and 
Astrophysics, CA 95064, USA 

Universitat Rostock, Institut fur Physik, D-18051 Rostock, Ger- 
many 



not Ipaumon and Guillotll2004l) , and if heavy elements 
in the interior are distributed homogeneously. These 
are fundamental interior structure properties. If we 
are to describe the interior of Jupiter and other giant 
planets, we shape our imagination of interior structure 
into terms such as core mass and envelope metallic- 
ity and apply the simplest possible structure type that 
can quantitively answer our questions. For instance, 
the assumption of a completely homogeneous interior 
is sufficient to explain the observed mass and radius of 
a extrasolar giant planet, and an occasional additional 
constraint can be satisfied by assuming the presence of a 
core within a two-layer approach. For solar giant plan- 
ets however, a two-layer approach is not consistent with 
available observational constraints whe n non-empirical 



equati ons of state are applied (e.g. : see IChabrier et al 



1 19921) for Jupiter and Saturn, and lHelled et all (|2010l ) 
for Uranus and Neptune) , and hence we here model the 
interior by assuming three layers. This article is not 
about the difficulties we run into by relying on this ap- 
proach. It is about the resulting core mass and metallic- 
ity of Jupiter fij l3.ll). where we repeat some of th e ana- 
lyis published in Fortnev and Nettelmann ( 2010l ) (here- 
after FN10). We then continue in § 13.21 to investigate 
the specific open question if the assumption of an ho- 
mogeneous heavy element dstribution can be subject 
to free choice in the case of recently published ab-initio 
da ta based Jupiter mod e ls, as it has been suggested 



Militzer and Hubbard! |2009l) . In (§15) we show the 



degeneracy of solutions raised by ignoring the octupole 
moment of the gravity field, J4 equivalent to knowing 
the Love number fc 2 of an extrasolar planet, but not 
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2 Modeling giant planets in isolation 

2.1 Observational constraints 

Interior models should reproduce the observed param- 
eters. These are, first of all, mass M p and radius R p 
of the planet. For Jupiter, Saturn, and Uranus, the 
radius is the observable equatorial radius R eq taken 
at the 1 bar pressure level. Below that the interior 
is adiabatic. For close-in extrasolar giant planets, the 
observed transit radius occurs at 4—6 orders of mag- 
nitude lower pressures than the onset of the adiabatic 
region, and a proper outer boundary condition has to be 
provided by a model atmosphere. Calculations for par- 
ticular giant planets predict the radiative/convective 
layer boundary deep in the atmosph ere between 100 



and 1000 bars ( Showman et al. 20081) . The mean he 



hum abundance should equal that of the protostellar 
cloud where the planet formed from. However, the at- 
mospheres of Jupiter and Saturn are depleted in he- 
lium compared to the solar valu e Y = 0.270 ± 0.005 
( Bahcall and Pmsonneault]ll995| ). which motivates the 
introduction of a He-poor outer envelope and a He-rich 
inner envelope. The metallicity Z in the atmospheres 
of the solar planets is enh anced above the solar value 
Z Q — 0.015 ( Loddersl 20031) as derived from abundance 
measurements of single species (e.g.; C, N, S, Ar, Xe). 
Solar planets are fast rotators, which gives rise to a 
gravity field deformation as expressed by the gravita- 
tional moments J2, J a, and Jg. Equivalent constraints 
for extrasolar planets are tidal Love numbers km- Pa- 
rameters related to thermal or orbital evolution can give 
additional constraints but are not addressed in this ar- 
ticle. 

2.2 The Three-layer structure assumption 

For Jupiter and Saturn, the simplest structure type 
that is consistent with the constraints described above 
has three layers: a core of rocks and/or ices, and two 
envelopes that are convective, adiabatic, and homoge- 
neous, but differ in the mass fraction of helium (Y) and 
metals (Z). We choose the m etallicities Z\ and j?J]] to 



adjus t J 2 and J4. Of ten (e.g; Guillot 1999 ; lHori et al 



20081) but not always (jMilitzer et al.ll2008D models that 
reproduce Jupiter's J2 and J4 are found to reproduce 
also Jq within the obcrvational error bar. Physical rea- 
sons for a discontinuity in He or metals can be phase 
transitions, H/Hc phase separation, and in particular 



1 Index 1 refers to the outer envelope, and index 2 to the inner 
envelope. 



for metals the process of planet formation, but its lo- 
cation is essentially unconstrained, and hence the tran- 
sition pressure P1-2 between the envelopes a free pa- 
rameter. In the absence of these particular physical 
processes, gaseous planets are likely to have a homoge- 
nous envelope because of convection. If h owever convec- 



tion i s inhibited as suspected for Uranus (jPodolak et al 



19911) a series of double-diffusive layers can develop cre- 
ating a large-scale compositional gradient. Such com- 
plicated models are not considered here. 

2.3 Quasi-adiabat and equation of state 

At layer boundaries, pressure P and temperature T are 
required to transit continously. The P — T — p pro- 
file, where p is mass density, along the quasi-adiabat 
depends on the He abundances and metallicities cho- 
sen, and on the equations of state (EOS) of the un- 
derlying materials allowed for. If the envelopes differ 
in He mass fraction (Y\ < Y2) or metallicity [Z\ 7^ 
Z2), then entropy and density transit discontinously, 
so we call the internal P — T profile quasi- aAiahtdic. 
At a layer boundary, a boundary layer can develop 
with change in temperature across or conductive heat 
transport, causing a warmer interior at higher entropy. 
The effect of a boundary lay er has been included by 
Fortnev and Hubbard! ( 2003 ) for the inhomogeneous 
evolution of Saturn; but not yet in structure calcula- 
tions. Since the gravity data probe the internal P — p 
relation, the expected influence of a warmer interior on 
the structure would be a higher deep envelope metal- 
licity. 

It is convenient to represent elements heavier than 
He (metals) by an equation of state of water (e.g. ; 
Fortnev and Hubbardll2003HSaumon and Guillotll2004h 
as als o supported by planet formation theory (jHelled et al.[ 
20081), and to assume the core be made of rocks 



2004; 



Nettelmann et al 



Saumon and Guillot 



e.g.: iFortnev and Hubbard! l2003t ISaumon and Guillot 



200 



4). 



20081 ) or of water ice (e.g.; 
Jupiter-sized giant plan- 



ets are giant in size because they are predominantly 
composed of the light e lements H and He. For de - 
tails of the EOS, see ISaumon and Guillot (2004); 



Nettelmann et al 
( 20081) . 



(|2008l) . FN10, and iMilitzer et al 



2.4 Calculation of core mass and metallicity 

In order to obtain core mass and metallicity we inte- 
grate the 1-dimensional equations of mass conservation, 
dm/dl = 4irl 2 p(l), and of hydrostatic equilibrium, 



1 dP _dU 
plj)~dT ~~dl 



(1) 
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along the quasi-adiabat for given outer boundary condi- 
tions and abundances of helium and metals. The coor- 
dinate I parametrizes surfaces of constant total poten- 
tial U. In spherical symmetry, i.e. in the absence of ro- 
tation, I equals the radial coordinate r. The most often 
theory used to calculate the gravity field deformation 
ac cording to shape deforma t ion is the Theory of Figures 
by Zharkov and Trubitsvnl ( 1978 ) , but alternativ e the- 
ories can be applied as well (e.g.; lHubbardlll984l ). The 
gravity field deformation on or exterior to the surface 
(p(P = lbar) w 0) can be described by multipole ex- 
pansion into spherical harmonics, where the expansion 
coefficients are the gravitational moments 



variety of currently competing hydrogen equations of 
state. The total mass of metals M core + Mz is better 
constrained and amounts to 14— 38 Mq (more precisely, 
in the order of the above EOS listing: 14 - 24, 20 - 24 
and 28 - 32, 17 - 37, 33, 38, and 35 M e ). Thus to 
our current knowledge 4.4 — 12% of Jupiter's mass are 
heavy elements, which is an enrichment factor of 3 — 8 
over solar mctallicity. 



J2n 



1 



M p R% 



d i rp(r,6)r' tn P 2n (cos6) 



(2) 



They are intergrals over the internal mass density 
within the volume V of the planet. Therefore, we can 
use the metallicities Z\ and Z 2 , which affect p(l) along 
the quasi-adiabat, to adjust the moments J 2 and J4. In 
general Z\ 7^ Z 2 . The core mass M core is then used to 
satisfy the inner boundary condition m(l = 0) = for 
given Mp(Rp). 



3 Jupiter 

3.1 Predicted core mass and envelope mass of metals 

Figure [T] is a modified version of Fig. 1 in FN10 and 
shows model results for the core mass and envelope 
mass of metals of Jupiter predicted by different equa- 
tions of state. These are DFT-Mrj§ |Militzer et al 
20081 ). LM-REOS based on FT-DFT-ME0 with differ- 
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Fig. 1 Model results for the core mass (M core ) and mass 
of metals (Mz) in the envelope(s) of Jupiter. Each poly- 
gon embraces the solutions found using a particular set of 
equations of state as labeled, for the components H, He, 
and metals. Models using SCvH-i (SG04) and Sesame- 
p as sume Z\ = Z2 and Y\ < Y < Y\ and are taken 
from lSaumon and Guillot) (I2004T ), SCvH-i (3L) based mod- 
els have Z\ ^ Z->. |Guillotlll999h . the Sesame-K04 model has 



Zi < Z 2 (lKerlev(l2004h as have LM-R EOS based models 

and the DFT-MD 
Y = 0.24 



ent EO S representing metals, n amely He4 and H2O 
REOS dNettelmann et al-lfiooih. SCvH-i and Sesame 



JN 

p (|Guillotlll999HSaumon and Guillotll2004l). Sesame-92 



(updated from iNettelmann et al.ll200l 
labeled model is a two-layer model with Y\ 
throughout the envelope. Adjusting this model for the pro- 
tosolar mean value Y = 0.27 by replacing 9 M® of metals 
by He yields a metal free envelope as indicated. This fig- 
ure shows < Mcore < 18 M® , < M z < 38 M© , and 

"U4 < Mcore + Mz < 38, meaning that Jupiter's interior 
(the original Sesame-EOS for H and He ( Lyon a nd JohnsonF — . , ,.7 

- 1 1 „ „ , ; structure is not well known today and that models depend 

1223) together with H 2 0-REOS), and Sesame-K04 (a on the equation of gtate preferr J 

Since core mass and envelope mass of metals are de- 
termined by the requirement to find internal density 
distributions that reproduce M p (R p ), J 2 and J4, the 
results reflect the compressibility of the H/He subsys- 
tem along the quasi-adiabat. It is then particularly 
strange that similar methods used to calculate H and 
He EOS yield contrary results on M core and Mz, as 
it is the case for DFT-MD and LM-REOS based on 
FT-DFT-MD 



revised H and He Sesame-EOS wi th minor cons tituent 
EOSs from the Panda-Code; see iKerlevI koo4 ) . The 
Jupiter models are published in the same references 
apart from t he Sesame-9 2 base d models, which were 
calculated in Nettelmann (2009). For each particular 
equation of state, numerous solutions are found due to 
uncertainties in the constraints applied. Obviously, a 
large uncertainty in Jupiter's core mass (0 — 18M®) and 
envelope metallicity (0 — 38 M ffi ) arises from the broad 

2 Density Functional Theory for the electrons and Molecular Dy- 
namics simulations for the ions in the warm dense matter regime 

3 In addition to DFT-MD, inclusion of finite temperature effects 
on the electronic subsystem, which affect the forces onto the ions 



sec 



Nettelmann et al. 2008 



2008HKietzmann et al.ll2007t iFrench et al 



Hoist et al 



2009, for de- 



tails). 
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3.2 The freedom to assume a discontinuity in heavy 
elements 



In § 13.11 some of the models were based on the assump- 
tion Z\ = Zi, while others (i.e. those using SCvH-i-99, 
Sesame-K03, Sesame-92, LM-REOS) were not. Here 
we show that this assumption is not for every EOS 
subject to free choice, if the model is required to meet 
J2, J4, and Jq. In particular, LM-REOS based Jupiter 
models require Z\ «C Z2, while DFT-MD based models 
Z\ w Z2. 

The convergent procedure used for LM-REOS based 
models is illustrated in Figj2j Along the thick solid line, 
Z\ = Z-2- On that line, there is one model that matches 
J2. However, Ji/10~ 4 = —6.04 for this model is 8a afar 
from the observed J4 value (see figure caption for de- 
tails). Hence the convergent procedure follows the thin 
black line in the direction of the arrows until both J2 
and J4 are met. Along that path, the core mass de- 
creases from 12 Mq down to 1.9 M®. The converged 
solution has Z\ < Z2. This behavior applies to every 
LM-REOS based Jupiter model, and hence the choice 
Z\ = Zi cannot be made unless the interior does not 
rotate rigidy but inst ead the observed J4 v alue is influ- 
enced by deep winds (jMilitzer et al.ll2008l ). 

On the other hand, as Fig. [2] shows, Z\ <C Z2 is 
not a favorite option for the DFT-MD based Jupiter 
model, since that would imply Z\ <C Zq or even 
zero. This would contradict heavy element abundance 
measurements in Jupiter's atmosphere, which indicate 
Z\ > 2 x Zq. Therefore, Z\ « Z2 is not a free choice 
for the DFT-MD based model. Consequently, the very 
different Jupiter core masses obtained using these two 
ab-initio EOS are not primarily due to different assump- 
tions about the distribution of heavy elements. In order 
to resolve this problem, a comparison of the H/He adi- 
abats is highly desirable. Furthermore, isentropic com- 
pression experiments in the 0.5 — 4 Mbar regime where 
J2 and J4 are most sensitive to mctallicty would greatly 
help to discriminate between competing Jupiter mod- 
els. For the other EOS considered here, both the as- 
sumptions Z\ = Z2 or Z\ 7^ Z2 give acceptable Jupiter 
models and hence are a matter of free choice. 




Fig. 2 Converged and intermediate envelope metallici- 
ties Z\, Z2 of a three-layer Jupiter model based on LM- 
REOS. The converged model matches the most recent valu e 
J4/IO -4 value (grey numbers) of —5.87 (|jacobsonl 120031 ') . 
has M core = 1.9 (black numbers) and Z\ = 1.6%. 
This model was obtained by an interative procedure that 
produced intermediate solutions along the solid black line. 
Starting from Z\ — Z2 = 1%, the envelope metallic- 
ity is increased along the Zi — Zi line of homogeneous 
metallicity (thick black solid) until J2 is matched. Such a 
model with LM-REOS would have M core = 12.0 M e and 
J4/10~ 4 = —6.04. Then iteratively J4 and J2 are fitted by 
allowing for Z\ 7^ Z2 as indicated respectively by grey circles 
and black circles until outer and inner envelope metallicities 
are found that match both J2 and J4 (converged model). All 
models along the dashed black line match J2 , and all models 
along the dashed grey line match J 4. This procedure can- 
not be applied to bring the J4 value of the D FT-MD based 
Jupit er model (grey square with error bars) (jMilitzer et alj 
l2008h in agreement with the observed J4 value, since obvi- 
ously Z\ <C Zq would rapidly occur. 



planets, currently available gravity data are M p and R p 
only. A potentially observational constraint equivalent 
to J2 has been suggested to be the tidal Love num- 
ber &2 (jRagozzine and Wolfll2009l ). which measures the 
ability of a planet to develop an elliptic deformation in 
response to a tidal perturber such as the close-by parent 
star. In this section, we present implications of know- 
ing a precise &2 value of the Neptune-sized extrasolar 
planet GJ 436b. 



4 Extrasolar planet models 

4.1 J 2 and &2 

As we have seen in §[3l Jupiter models that use one par- 
ticular EOS can have various resulting {M core , Z±, Z2} 
triples. Without the constraint imposed by J4 how- 
ever, the variety of solutions would be immense, includ- 
ing M core > 18 M® and Mz > 38 M®. For extrasolar 



4.2 Core mass of GJ 436b models for given &2 

About 30 light years away from Earth, the Hot Neptune 
G J 43 6b (M p = 23.17 M e , R p = 4.22 R @ ) ((Torres et al ' 



20081) orbits the M-star GJ436. Placed in a mass-radius 
diagram, this planet is located close to theoretical M-R 
relations of warm water planets. Accordingly, interior 
structure models assuming an iron-silicate core, a water 
layer, and a H/He envelope allow for a composition of 
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95% water and 5% H/He jFigueira et al.l2009f ). but also 
for a water- less composition. H/He envelope models, 
whether dry or not, are found to have 0.02 < k 2 < 0.2 



( Nettelmann et al. 2010t) . We here assume k 2 = 0.2 



and search for the core mass and water content that 
satisfy this additional constraint. We allow water be 
mixed homogeneously into the H/He envelope (two- 
layer model) or confined to a deep water layer (three- 
layer models) and combinations in between (Neptune- 
like three-layer models). Figure [3] shows the result 
in comparison with LM-REOS based Jupiter models 
forced to meet the observed J 2 value, but not J4. 



0.1 



Jupiter models 

J 2 =0.014697 



J4/IO =-6.05 




envelope M z / M p 



0.5 

envelope M z / M p 



Fig. 3 Jupiter models matching the measured J2 value, 
and GJ 436b models matching an arbitrary value. Max- 
imum core mass (no water layer) and minimum core mass 
(no rock core) solutions are found. Accurate (~ 1%) knowl- 
edge of Jupiter's J4 reduces the uncertainties in M core and 
Mz {filled circle solution), while GJ436b's core mass is un- 
determined if &2 = 0.2 and k± remains unknown. 



5 Conclusions 

Standard Jupiter interior models predict a core mass of 
— 18 M®, an envelope mass of metals of — 38 M®, 
and a total mass of metals of 14 — 38 Mm . These uncer- 



tainties have increased over the past ten years (|Guillot 
1999). when they were thought to arise from the pos- 
sible presence of a plasma phase transition in hydro- 
gen. From the behavior of model solutions discussed in 
this article, we conclude an unsufficient understanding 
of the H/Hc EOS in the 0.5—4 Mbar pressure regime, 
where J 2 and J4 are most sensitive. Equivalent to 
Jupiter models, measuring the k 2 value of an extrasolar 
planet would reduce possible solutions to a linear rela- 
tion between core mass and envelope mass of metals, 
as shown for GJ 436b. 

The author kindly thanks R. Redmer, J.J. Fortncy, 
T. Guillot, and D.J. Stevenson for insightful discus- 
sions, U. Kramm for providing the k 2 values, and the 
anonymous referee for valuable comments. 



In both cases, the solutions are located along a 
straight line (solutions for different k 2 or J 2 values 
would span parallel lines). For both planets, and hence 
in general, the upper limit in core mass is obtained for 
homogeneous envelope models (compare § I3.2[) . and the 
upper limit in envelope mass of metals (here water) by 
zero-mass core models. In between there are three-layer 
models with various mctallicitics in the two envelopes. 
As indicated for Jupiter, solutions below the filled circle 
are also in agreement with J4. 

We conclude that the constraint imposed by k 2 or 
J 2 reduces the degeneracy of models, but further con- 
straints are needed to determine the core mass of a 
planet within a few percent of total mass (Jupiter: 
6%). Cooling curve calculations might help to constrain 
the composition, since they depend on specific heat 
which significantly varies between relevant bulk materi- 
als. Unfortunately, evolution calculations are subject to 



uncertainties such as tidal interactions (Leconte et al 



2010) and stella r irradiation and th e response of the 
atmosphere (e.g. Valencia et al. 2010[) . 
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